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1.  The  dependence  of  observed  surface  wave  amplitudes 
on  explosion  yield  and  the  characteristics  of  the  emplacement 
medium  is  reviewed.  Surface  wave  data  compiled  by  others  is 
used  in  conjunction  with  a new  data  set  consisting  of  Airy 
phase  amplitudes  from  the  WWSSN  stations  ALQ  and  TUC. 

2.  The  theoretical  dependence  of  surface  wave  amplitude 
on  the  important  controlling  parameters  is  discussed.  Atten- 
tion is  directed  to  the  assumptions  of  the  theoretical  models 
and  their  consequences. 

3.  Preliminary  results  are  presented  for  the  develop- 
ment and  testing  of  a deconvolution  technique  for  analyzing 
short  period  teleseismic  recordings  of  explosions. 
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I . SUMMARY 


The  objective  of  our  research  program  is  to  examine 
the  parameters  that  affect  the  seismic  signals  from  under- 
ground explosions.  Our  attention  is  primarily  directed  to 
those  features  of  the  seismic  waveforms  that  reliably  indicate 
the  explosion  yield.  Our  research  program  includes  empirical 
studies  of  the  available  data,  experimental  studies  using 
small  charges  to  simulate  explosions  and  the  development  and 
application  of  theoretical  and  numerical  methods.  Emphasis 
is  on  the  latter.  In  particular,  we  are  applying  techniques 
for  numerically  simulating  the  far-field  signals  from  both 
contained  and  cratering  underground  explosions.  The  numeri- 
cal simulation  techniques  represent  a synthesis  of  the  finite 
difference  methods  for  computing  ground  motion  in  the  near- 
source large  displacement  regime  and  the  efficient  wave  propa- 
gation techniques  of  theoretical  seismology. 

During  the  third  three-month  period  of  our  present 
contract,  our  research  has  been  conducted  in  a number  of 
areas.  These  areas  and  the  major  results  obtained  in  each 
are  summarized  in  the  following  paragraphs. 

A considerable  amount  of  effort  has  been  directed  to 
the  study  of  explosion  generated  surface  waves  during  this 
quarter.  Our  objective  is  to  understand  the  dependence  of 
surface  wave  amplitude  on  explosion  yield,  the  characteris- 
tics of  the  emplacement  medium  and  the  characteristics  of  the 
source-receiver  travel  path.  The  most  commonly  used  measure 

of  surface  wave  amplitude  is  M_  and  we  have  studied  the  M 

s s 

data  compiled  by  others.  However,  we  have  also  developed  our 
own  data  base  for  the  Nevada  Test  Site  (NTS)  using  Airy  phase 
measurements  from  the  WWSSN  stations  ALQ  and  TUC.  These 
stations  are  within  1000  km  of  NTS  and  have  the  advantage  of 
recording  events  over  a large  yield  range.  For  events  in 
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common,  we  find  the  relative  amplitudes  from  ALQ  and  TUC  to 

be  consistent  with  the  best  M„  data. 

s 

Comparing  events  in  different  local  materials  at  NTS, 
we  find  that  the  M coupling  is  about  0.1  to  0.3  units  higher 
for  events  below  the  water  table  at  Pahute  Mesa  than  for 
comparable  everts  at  Yucca  Flat.  While  there  are  only  a few 
events  above  the  water  table  at  Pahute  Mesa,  they  clearly 
couple  more  weakly  into  Mg  than  those  in  the  saturated  mate- 
rials. The  lowest  coupling  events  in  our  population  are  those 
in  dry  tuffs  at  Yucca  Flat.  While  the  scatter  is  considerable, 
these  events  are  0.5  to  0.6  units  lower,  on  the  average,  than 
saturated  events  at  Yucca  Flat.  The  saturated  tuff  events  in 
the  tunnel  beds  at  Rainier  Mesa  seem  to  couple  about  the  same 
as  the  saturated  tuff  events  at  Yucca  Flat.  Finally,  the 
PILEDRIVER  explosion  in  granite  seems  to  couple  a small  amount 
higher  than  events  in  the  highest  coupling  population,  the 
Pahute  Mesa  explosions  below  the  water  table. 


Theoretical  studies  of  surface  wave  generation  and 
propagation  are  of  considerable  value  for  organizing  and 
understanding  the  empirical  data.  It  is  important,  therefore, 
to  carefully  set  forth  the  theoretical  results,  making  clear 
the  assumptions  being  made  and  organizing  these  results  in 
the  proper  context  for  comparing  with  data.  As  most  commonly 
used  for  explosion  studies,  the  theory  assumes  a spherically 
symmetric  point  source  representation  for  the  source  and  a 
plane-layered  model  for  the  path.  This  theory  is  fine  as 
far  as  it  goes,  but  as  we  study  surface  waves  more  closely, 
certain  limitations  become  apparent.  Among  the  most  important 
of  these  is  the  need  to  account  for  the  fact  that  the  explo- 
sion source  is  certainly  more  complex  than  a spherically 
symmetric  point  source.  It  is  also  important  to  account  for 
the  fact  that  the  crustal  structure  in  the  source  vicinity 
where  the  surface  waves  are  originally  excited  is  often 
rather  different  than  that  for  the  rest  of  the  path  which 
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controls  the  dispersion.  Making  a few  reasonable  assumptions, 
our  theoretical  formulation  can  be  adjusted  to  account  for 
these  limitations. 

The  surface  wave  amplitude-yield  relationship  is  con- 
trolled by  several  interlocking  effects.  Explosion  yield 
changes  are  almost  always  accompanied  by  changes  in  the  burial 
depth  and  the  properties  of  the  emplacement  material.  In 
turn,  the  burial  depth  has  two  separate  effects;  one  due  to 
the  changing  overburden  pressure  influence  on  the  source  func- 
tion and  one  due  to  the  depth  effect  on  surface  wave  excita- 
tion. If  the  latter  effect  is  isolated,  we  find  that  surface 
wave  excitation  is  a slightly  increasing  function  of  depth. 

That  is,  if  all  other  factors  were  fixed,  M -log  yield  curves 
would  have  a slope  of  about  1.06  for  NTS  tuff  and  slightly 
greater  than  unity  for  NTS  granite.  Of  course,  all  other 
factors  are  not  equal,  but  these  too  can  be  isolated  to  deter- 
mine their  relative  influence. 

If  we  assume  a spherically  symmetric  point  source 
representation  for  the  explosion,  we  find  surface  wave  ampli- 
tude to  be  strongly  dependent  on  the  local  shear  modulus. 

In  fact,  we  find  M ss  log  [y^*  ¥ ] , where  is  the  (equivalent 

O 

elastic)  source  which  is  found  to  be  nearly  independent  of  . 
period  for  periods  near  those  at  which  Mg  is  measured.  But 
we  know  the  source  is  not  spherically  symmetric.  In  fact,  it 
is  easy  to  argue  that  the  spherical  source  for  waves  traveling 
upward  from  the  source  is  different  from  that  for  waves 
traveling  downward.  How  different  we  do  not  yet  know,  though 
the  question  is  under  study. 

In  studying  the  relative  contribution  of  the  up  and 
downgoing  portions  of  the  wave  field  to  teleseismic  surface 
waves,  it  is  interesting  to  contrast  two  types  of  source 
regions.  In  one  the  source  is  in  a soft  material  like  tuff 
which  overlays  the  hard  basement  materials.  In  the  second 
the  harder  materials  (granites)  extend  nearly  to  the  surface 
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and  this  soft  rock-hard  rock  contrast  below  the  source  is 
absent.  For  the  latter  case  the  surface  waves  from  the  up 
and  downgoing  portions  of  the  source  are  in  phase  and  add. 
However,  when  the  sharp  contrast  is  present,  the  up  and 
downgoing  generated  surface  waves  are  out  of  phase  and  tend 
to  cancel.  One  interesting  result  is  then  that  if  the  upgoing 
waves  were  entirely  suppressed  or  scattered,  the  dependence 
of  teleseismic  surface  wave  amplitude  on  y is  much  reduced. 

In  addition  to  the  work  on  surface  waves  that  is  briefly 
summarized  in  the  previous  paragraphs,  we  have  been  working 
with  deconvolution  techniques  for  analyzing  short  period  body 
wave  recordings  of  explosions.  Teleseismic  body  wave  re- 
cordings can  be  viewed  as  recordings  of  the  source  generated 
waves  modified  by  a series  of  transfer  functions.  There  are 
many  combinations  of  events  and/or  recordings  that  allow  isola- 
tion of  one  or  more  of  these  transfer  functions.  Once  they  are 
known,  they  can  be  deconvolved  from  recordings  of  future  events 
to  isolate  the  contribution  of  the  source. 

A signal  deconvolution  technique  based  on  "homomorphic 
deconvolution"  methods  described  in  the  literature  has  been 
programmed  and  is  being  tested  on  data.  The  success  of  the 
technique  must  really  be  judged  on  the  transfer  functions  it 
gives  and,  particularly,  on  how  stationary  they  are  from 
experiment- to-exper iment . 


II.  INTRODUCTION 


The  primary  objective  of  the  research  program  is  to 
systematically  examine  the  parameters  that  control  magnitude- 
yield  relationships  for  underground  nuclear  explosions. 

During  the  third  three-month  period  of  this  contract  research 
has  been  conducted  in  the  following  areas: 

1 . A review  of  the  dependence  of  observed  surface 
wave  amplitudes  on  explosion  yield  and  the 
characteristics  of  the  emplacement  material 
(summarized  in  Section  III) . 

2.  A theoretical  study  of  the  dependence  of 
surface  wave  amplitude  on  explosion  yield, 
the  characteristics  of  the  emplacement 
material  and  the  characteristics  of  the 
travel  path  (summarized  in  Section  IV) . 

3.  Development  and  testing  of  a deconvolution 
technique  for  analyzing  short  period  tele- 
seismic  recordings  of  explosions  (summarized 
in  Section  V) . 

The  remainder  of  the  report  is  devoted  to  technical 
discussions  of  results  in  these  three  areas. 
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III.  SURFACE  WAVE  OBSERVATIONS 

3.1  INTRODUCTION 

We  wish  to  understand  the  dependence  of  surface  wave 
amplitude  on  explosion  yield  and  the  characteristics  of  the 
emplacement  medium.  The  surface  wave  amplitude  data  we  are 
using  come  from  two  sources.  First,  there  are  the  M data 
given  by  Eisenhauer  11976] . Second,  we  have  Airy  phase  or 
maximum  amplitude  measurements  from  two  WWSSN  stations  at 
less  than  1000  km  from  NTS.  In  this  section  we  first  sum- 
marize these  data  and  then  compare  them. 

The  reason  for  using  the  Airy  phase  data  is  to  attempt 
to  extend  the  surface  wave  data  base  to  as  many  different 
testing  areas  as  possible.  We  would  like  to  be  able  to  esti- 
mate the  relative  surface  wave  coupling  of  events  above  the 
water  table  and  at  Rainier  Mesa  compared  to  the  Pahute  Mesa 
and  Yucca  Flat  events  below  the  water  table.  There  are 
insufficient  Mg  data  for  thir,  purpose. 


3.2  Ms 


The  most  reliable  data  for  Mg  seem  to  be  those  compiled 

by  Eisenhauer  [1976].  The  data  of  interest  for  our  purposes 

are  those  for  NTS  explosions  at  or  below  the  water  table  or 

in  granite.  Only  a single  even.  (PILEDRIVER)  falls  in  the 

latter  category.  The  other  events  are  separated  according  to 

whether  they  were  at  Pahute  Mesa  or  Yucca  Flat.  Linear  least 

squares  equations  of  the  form  M ■ a + b log  W were  computed 

8 

for  each  of  the  two  populations.  The  results  are: 

Pahute  Mesa:  M ■ 1.07  log  W + 1.88, 

8 (1) 
Yucca  Flat:  M#  ■ 1.08  log  W + 1.70. 
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Since  yields  near  150  kt  are  of  primary  interest, 
intercepts  for  an  equation  of  the  form  M ■ a + log  W were 
derived  for  each  test  area  using  explosions  between  60  and 
415  kt.  There  were  14  events  in  the  reduced  Pahute  Mesa 
population  and  21  explosions  at  Yucca  Flat  in  this  yield 
range.  Assuming  this  unit  slope  relationship,  values  of 
the  intercept  (a)  can  be  computed  for  each  event.  If  we 
average  these  values  we  find  the  following  (a  is  the  standard 
deviation) : 

Pahute  Mesa:  Intercept  = 2.00,  o = 0.14, 

Yucca  Flat:  Intercept  = 1.87,  a = 0.13,  (2) 

PILEDRIVER:  Intercept  = 2.09. 

Another  useful  way  to  present  the  "average"  M for  each  testing 

s 

area  is  by  computing  the  value  for  a 150  kt  event. 


Equation  (1) 

Equation  (2) 

Pahute  Mesa: 

4.21 

4.18  + 0.14 

Yucca  Flat: 

4.05 

4.05  + 0.13 

PILEDRIVER: 

4.26 

4.26 

3.3  AIRY  PHASE  AMPLITUDE  FROM  WWSSN  STATIONS  ALQ  AND  TUC 

Data  from  the  WWSSN  stations  ALQ  and  TUC  were  collected 
for  a large  sampling  of  NTS  explosions.  Some  results  were 
reported  by  Bache,  et  al . , [1975a].  The  station  distances 
4 and  azimuths  are  summarized  below: 


L/ 
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ALQ 


TUC 


Test  Area 

Distance 

(km) 

Azimuth 

Distance 

(km) 

Azimuth 

Yucca  Flat 

900 

103° 

720 

136° 

Rainier  Mesa 

910 

103° 

730 

136° 

Pahute  Mesa 

930 

103° 

755 

135° 

PILEDRIVER 

900 

104° 

730 

137° 

These  values  are  for  events  near  the  center  of  each  test  area. 
Typical  seismograms  for  these  two  stations  are  shown  in 
Figure  1. 

The  peak-to-peak  amplitudes  of  the  maximum  cycle  in 
the  Rayleigh  wave  train  and  the  period  of  that  cycle  were  mea- 
sured for  each  event.  All  measurements  were  from  the  vertical 
component.  Since  the  propagation  path  lies  entirely  in  a con- 
tinental structure  and  the  instrument  passband  is  7 to  30 
seconds,  the  maximum  amplitude  is  associated  with  the  Airy 
phase.  The  period  of  this  phase  was  11.5  + 0.5  seconds  at 
ALQ  and  8+0.5  seconds  at  TUC.  A striking  characteristic  of 
the  data  was  the  similarity  of  the  waveforms  from  event-to- 
event  at  each  station.  Thus,  the  amplitude  measurements  con- 
sistently sample  the  same  portion  of  the  waveform. 

The  Airy  phase  amplitudes  were  converted  to  M values 

s 

by  using  the  formulas: 

ALQ:  = log  A + 2.72, 

(3) 

TUC:  - log  A + 2.17. 


The  constant  factors  in  these  equations  were  chosen  to  make 
and  Mg  be  about  the  same,  on  the  average,  as  Eisenhauer's 
Ms.  No  other  meaning  is  attributed  to  them. 


The  ALQ  and  TUC  Mg  values  may  be  analyzed  in  about  the 
same  way  as  were  the  M data  in  Section  2.2.  First,  the  two 
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Figure  1.  Long  period  vertical  component  recordings  of  three 
NTS  explosions  at  two  WWSSN  stations. 
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measurements  are  averaged  to  obtain  our  best  estimate  of  Mg 
from  these  Airy  phase  data.  We  then  separate  the  data  by 
test  area  and  compute  a linear  least  square  fit  to  each 
population.  The  regression  is  done  assuming  the  yield  values 
are  exact  and  the  standard  deviation  of  the  Mg  values  from 
the  mean  is  determined.  The  linear  best-fit  equations  for 
each  population  together  with  the  standard  deviations  are 
as  follows  (the  number  of  events  is  listed  in  parentheses) : 

Pahute  Mesa  below  the  water  table  (9) : 

M = 0.86  log  W + B,  a * 0.09, 

S 

Yucca  Flat  below  the  water  table  (30) : 

M = 1.17  log  W + B - 0.93,  a « 0.19, 

S 

Pahute  Mesa  above  the  water  table  (3):  ^ 

M « 1.21  log  W + B - 1.09,  o - 0.05, 

Yucca  Flat  above  the  water  table  (15) : 

M - 0.78  log  If  + B - 0.82,  o - 0.27, 
s 

Rainier  Mesa,  tunnel  shots  only  (10) : 

Mb  - 0.77  log  W ♦ B - 0.41,  o - 0.19. 

The  Mg  data  together  with  the  least  square  fitting  lines 
are  plotted  in  Figure  2.  As  expected  from  the  large  a values 
in  (4),  there  is  considerable  scatter  in  the  data,  especially 
for  the  Yucca  Flat  events. 

How  well  do  the  Airy  phase  data  correlate  with  the 
Eisenhauer  M.  estimates?  For  saturated  events  at  Pahute  Mesa 
there  are  seven  events  for  which  both  data  are  available.  For 
Yucca  Flat  there  are  eleven  common  events.  In  Figure  3 the 
Airy  phase  M estimates  from  Figure  2 are  plotted  versus 
Eisenhauer '8  Mg  for  the  common  events.  While  there  is  some 
> scatter,  we  see  that  the  two  M estimates  are  consistent. 

S 

Taking  the  residuals  (M  [Eisenhauer]  - M [ALQ  + TUC])  and 
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PAHUTE  MESA  BELOW  THE  WATER  TABLE 
YUCCA  FLAT  BELOW  THE  WATER  TABLE 


Figure  2.  The  M values  from  the  WWSSN  stations  ALQ  and  TUC  are  plotted  versus 
exploiion  yield.  The  explosions  are  divided  into  different  popula- 
tions according  to  the  gross  properties  of  the  emplacement  media. 
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averaging,  we  find  the  mean  to  be  0.04  with  a standard 
deviation  of  0.10  Mg  units.  For  the  Pahute  Mesa  events 
taken  separately,  the  mean  residual  is  0.01  and  it  is  0.06 
for  the  Yucca  Flat  events.  The  small  mean  is  really  a test 
of  our  selection  of  the  constants  in  (3) . The  important 
value  is  the  o of  0.10. 

E[  [ 

3.4  THE  RELATIVE  M COUPLING  OF  EVENTS  IN  DIFFERENT 

MATERIALS  AT  NTS 

I 

We  are  interested  in  the  dependence  of  surface  wave 
amplitude  on  explosion  yield  and  on  the  characteristics  of 
the  emplacement  medium.  We  take  Mg  to  be  the  measure  of 
surface  wave  amplitude.  It  is  important  to  recognize  that 
it  is  generally  impossible  to  entirely  separate  yield  effects 
from  those  due  to  burial  depth  and  local  material  property 
variations.  While  Eisenhauer's  Pahute  Mesa  and  Yucca  Flat 
events  below  the  water  table  seem  to  be  populations  for 
which  the  local  material  property  variations  are  minimized, 
burial  depth  is  directly  proportional  to  yield  for  most 
events.  Therefore,  the  Mg  - log  W curves  in  (1)  also  incor- 

_ porate  M - depth  effects. 

s 

The  Airy  phase  data  discussed  in  Section  3.2  were 
collected  to  attempt  to  estimate  the  relative  surface  wave 
coupling  for  different  emplacement  materials.  Once  again, 
it  is  impossible  to  isolate  the  coupling  effects  since  the 
yield  and  depth  range  are  generally  different  for  different 
source  materials.  It  is  necessary  to  keep  this  in  mind  when 
comparing  Mg  values  from  the  different  test  areas. 

The  M data  in  (1)  show  M - log  W to  have  a slope 
s s 

of  about  1.07  for  events  in  saturated  tuffs  and  rhyolites. 
Considering  the  scatter,  the  data  of  (4)  and  Figure  2 do 
I not  contradict  the  conclusion  that  the  Mg  - log  Y slope  is 
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one  or  a bit  larger.  Further,  the  results  of  Figure  3 give 
confidence  that  the  Airy  phase  amplitudes  can  be  used  to 
estimate  the  relative  M_  coupling  for  explosions  in  different 
areas. 

Eishenhauer  computed  the  average  value  of  M - log  W 

D 

for  events  in  a limited  yield  range  (60  to  415  kt)  for  the 
Pahute  Mesa  and  Yucca  Flat  events  below  the  water  table.  The 
results  were  given  in  (2) . We  can  do  the  same  for  each  group- 
ing of  events  discussed  in  Section  2.3.  The  results  are  as 
follows : 

Pahute  Mesa  below  the  water  table  (9  events) , 

Intercept  = 2.06,  o = 0.11. 

Yucca  Flat  below  the  water  table  (21  events) , 

Intercept  * 1.78,  a = 0.17. 

Pahute  Mesa  above  the  water  table  (3  events) , 

Intercept  = 1.64,  o = 0.17. 

Yucca  Flat  above  the  water  table  (15  events) , 

Intercept  ■ 1.21,  a - 0.28. 

Rainier  Mesa  tunnel  explosions  (10  events) , 

Intercept  = 1.75,  o = 0.19. 

Yucca  Flat  Paleozoic  rock  (3  events) , 

Intercept  = 1.69,  a = 0.27 

PILEDRIVER: 

Intercept  = 1.95. 

These  results  together  with  those  from  (2)  are  presented 
graphically  in  Figure  4.  We  see  that  the  M and  Airy  phase 
data  agree  that  the  M coupling  is  about  0.1  to  0.3  units 
higher  for  the  events  below  the  water  table  at  Pahute  Mesa 
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than  comparable  events  at  Yucca  Flat.  There  are  only  a few 
events  above  the  water  table  at  Pahute  Mesa,  but  they  clearly 
couple  more  weakly  into  Mg  than  those  in  saturated  materials. 
The  lowest  M events  are  those  in  dry  tuffs  at  Yucca  Flat. 

The  scatter  is  quite  large  for  these  low  yield  events,  but 
M - log  W is  0.57  lower,  on  the  average,  than  for  the  satu- 

S 

rated  events  at  Yucca  Flat,  a number  much  larger  than  the 
standard  deviation  of  the  data.  The  saturated  tuff  explosions 
in  the  tunnel  beds  at  Rainier  Mesa  seem  to  couple  about  the 
same  as  the  saturated  tuff  events  at  Yucca  Flat.  The  same 
can  be  said  about  the  three  events  identified  as  being  in 
Paleozoic  rock  at  Yucca  Flat.  Since  these  events  were  deto- 
nated 23  meters  or  less  from  the  tuf f-Paleozoic  interface 
[Springer  and  Kinnaman,  1971],  the  identification  of  these 
events  as  being  in  paleozoic  rock  is  somewhat  ambiguous  and 
we  are  probably  not  seeing  the  true  differences  between  events 
in  saturated  tuff  and  the  paleozoic  rocks.  Finally,  the 
PILEDRIVER  event  seems  to  couple  like  the  highest  coupling 
population,  the  Pahute  Mesa  events  below  the  static  water 
table. 
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IV.  THEORETICAL  SURFACE  WAVE  GENERATION  BY  EXPLOSIONS 


INTRODUCTION 


In  this  section  we  first  outline  our  procedure  for 
computing  theoretical  seismograms  for  underground  explosions, 
making  clear  the  assumptions  being  made.  The  theoretical 
dependence  of  surface  wave  amplitude  on  explosion  yield, 
source  depth  and  properties  of  the  emplacement  material  is 


then  discussed.  We  also  address  the  dependence  of  Mg  on 


properties  of  the  travel  path. 


OUTLINE  OF  THE  THEORY 


The  computational  method  being  used  for  generating 
theoretical  surface  wave  seismograms  is  based  on  that  out- 
lined by  Harkrider  [1964,  1970],  The  formulation  of  the 
theory  is  entirely  in  terms  of  linear  elasticity.  The  effect 
of  anelastic  attenuation  is  included  via  an  empirically 
determined  Q operator.  The  explosion  source  is,  of  course, 
highly  nonlinear.  For  using  the  surface  wave  theory  it  is 
•then  necessary  to  represent  this  source  in  terms  of  an 
"equivalent  elastic  source."  For  spherically  symmetric 
explosions  this  equivalent  elastic  source  is  the  reduced 
displacement  potential. 


The  most  elementary  application  of  the  theory  is  to 
assume  that  the  entire  source-receiver  travel  path  can  be 
modeled  as  a single  multilayered  medium.  The  explosion 
source  is  assumed  to  be  represented  by  a reduced  displace- 
ment potential.  A great  deal  can  be  learned  about  explosion 
generated  surface  waves  from  this  theoretical  formulation, 
but  there  are  several  problems  to  be  addressed.  Among  these 


Say  we  were  interested  in  surface  waves  from 


NTS  to  a fixed  seismic  network.  The  local 
geology  changes  from  place  to  place  within 
NTS,  but  most  of  the  source-receiver  path 
remains  the  same.  If  we  are  restricted  to 
a single  path  model,  there  is  no  consistent 
way  to  deal  with  explosions  in  different 
local  media. 

We  know  that  the  explosion  is  not  a spheri- 
cally symmetric  source  of  elastic  waves 
embedded  in  one  layer  of  a layered  earth. 
The  question  is,  how  severe  is  the  effect 
of  assuming  that  it  is? 

The  previous  question  carries  with  it  a 
whole  set  of  subquestions  about  different 
aspects  of  the  behavior  in  the  source 
region.  Among  the  nonspherical  effects 
that  are  potentially  significant  for  sur- 
face wave  coupling  are: 

a.  The  nonlinear  region  intersects 
the  free  surface  for  most  events. 

The  theory  assumes  elastic  reflec- 
tion by  the  free  surface. 


example,  the  overburden  pressure 
changes  significantly  from  top  to 
bottom  of  this  region. 

In  the  remainder  of  this  section  we  first  outline  the  elemen- 
tary theory  and  then  the  modifications  employed  to  deal  with 
the  first  problem  listed  above.  The  second  and  third  problems 
are  addressed,  at  least  in  a tentative  way,  in  later  sections, 
particularly  Section  4.6. 

Assume  that  the  entire  travel  path  is  represented  by  a 
single  multilayered  earth  model  and  the  source  is  represented 
by  a reduced  displacement  potential.  Then  the  Fourier- 
transformed  vertical  component  of  the  fundamental  mode  Rayleigh 
wave  is 


w(r,u>) 


A 


iu)4iry  4*  (w) 
s 


Ks  ~R 


(5) 


where  V (w)  is  the  reduced  displacement  potential,  y is  the 
shear  modulus  of  the  source  layer,  r is  epicentral  distance 
and  cR  is  the  phase  velocity  of  the  Rayleigh  wave  component 
of  frequency.  The  AR  is  the  amplitude  response  of  the 
layered  medium  and  is  independent  of  the  source  depth.  The 
depth  excitation  factor  is  K_  which,  in  the  notation  of 

S 

Harkrider  [1964],  may  be  written 


uMh) 


K * 


w_ 


1 

U 


a*  (h) 

S VCR  ' 


(6) 


where  h is  the  source  depth,  u*  and  w*  are  radial  and  vertical 
velocities  and  is  the  normal  stress.  The  subscripts  o and 
s refer  to  the  free  surface  and  the  source  layer.  The  two 
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terms  in  (6)  are  tabulated  for  several  source  structures  by 
Harkrider  [1970],  To  account  for  anelastic  attenuation  and 
sphericity,  we  multiply  (5)  by 


(ae  sinA) 


where  y is  an  empirically  determined  attenuation  factor,  A 
is  range  in  degrees  and  ae  is  the  radius  of  the  earth. 

To  use  (5)  we  must  have  an  average  earth  model  for  the 
travel  path.  This  is  fine  for  the  dispersion  characteristics 
of  the  path  and  numerous  seismologists  have  demonstrated  that 
synthetic  seismograms  can  be  computed  that  look  very  much  like 
the  observed  seismograms  (from  which  data,  after  all,  the 
average  earth  models  are  constructed) . However,  a problem 
arises  when  the  material  in  the  local  source  region  is  not 
the  same  as  the  "average"  material  at  that  depth.  It  is  pos- 
sible to  ignore  this  problem  (and  many  people  have)  if  one  is 
dealing  with  a series  of  events  in,  more-or-less,  the  same 
source  material.  For  example,  one  could  construct  theoretical 
M -log  yield  curves  for  explosions  in  any  of  the  different  NTS 
source  materials  discussed  in  Section  II.  However,  the  problem 
is  impossible  to  ignore  if  one  is  attempting  to  compare  events 
in  close  proximity,  but  in  different  source  materials. 

This  problem  of  different  source  materials,  same 
average  path  is  one  that  has  interested  us  for  some  time. 

An  earlier  attempt  to  deal  with  the  difficulty  was  to  con- 
tinue to  use  a single  path  model,  but  to  let  the  source 
region  be  modeled  by  a spherical  plug  of  the  source  material. 

A modified  reduced  displacement  potential  incorporating  the 
effect  of  the  boundary  between  the  "source"  and  "global"  mate- 
rials is  then  used.  This  theory  was  worked  out  and  the  con- 
sequences were  explored  to  some  degree  by  Harkrider  and  Bache 
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and  the  results  were  reported  by  Bache,  et.  al . , [1975b]. 

While  the  results  are  interesting,  this  method  for  handling 
the  source  region  is  not  very  satisfactory.  The  method  we 
now  prefer  is  described  in  subsequent  paragraphs. 

A natural  way  to  approximate  the  different  source 
region,  common  average  path  situation  is  with  two  earth 
models,  one  for  the  source  region  and  one  for  the  remainder 
of  the  path.  For  this  we  must  be  able  to  account  for  propaga- 
tion of  Rayleigh  waves  across  the  boundary  between  the  two. 
Alewine  [1974]  gives  a correction  to  account  for  the  boundary 
that  is  based  on  the  work  of  McGarr  [1969]  and  McGarr  and  Alsop 
[1967] . Letting  the  source  and  receiver  portions  of  the  path 
be  denoted  by  subscripts  1 and  2,  the  equation  for  mixed  path 
Rayleigh  waves  according  to  these  authors  is: 


w = - iw4Try  $ 
s 


(-Ylrl-Y2r2) 


rl  + r2 


1/2 


ae  sin 


(A1  + A2) 


(8) 


The  transmission  coefficient,  T(u>),  is  derived  by 
assuming  the  total  horizontal  energy  flux  remains  constant 
during  the  transmission  of  the  Rayleigh  wave  across  the 
boundary.  The  coefficient  is  given  by 


where  U is  the  group  velocity  and  E is  the  total  potential  or 
kinetic  energy  normalized  to  uQ,  the  surface  vertical  dis- 
placement. In  using  this  coefficient  we  assume  that  there  is 
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boundary  is  far  enough  from  the  source  so  that  the  wave 
front  has  small  lateral  curvature.  Finally,  since  no  refrac- 
tion effects  are  included,  we  assume  that  the  wave  front  is 
parallel  to  the  boundary. 

How  do  we  use  (8)  to  study  NTS  explosions?  First, 
as  far  as  we  know,  the  only  difference  in  source  structure 
between  events  in  different  portions  of  NTS  is  in  the  top 
few  kilometers  of  the  structure.  We  need  the  freedom  to 
change  this  portion  from  granite  to  tuff  or  some  other  mate- 
rial. Then  for  the  long  periods  of  interest  for  M , the 

5 

boundary  is  nearly  transparent  and  T(w)  « 1.  Even  at  short 
periods  it  seems  reasonable  to  let  r^  ■+•  0 because  the  disper- 
sion and  attenuation  are  average  path  quantities  that,  to 
the  extent  they  are  derived  from  NTS  data,  already  incorporate 
any  mixed  path  effects  that  might  be  present. 

As  far  as  the  longer  period  surface  waves  and  Mg  are 
concerned,  the  practical  consequences  of  using  this  formula- 
tion are  simple  and  attractive.  The  depth  excitation  factor, 

K , is  computed  using  a structure  that  has  the  correct  mate- 
rial  at  the  source  depth.  There  is  then  no  ambiguity  in  the 

use  of  u . That  is,  the  material  for  7 is  the  same  as  the 
s 

material  in  the  multilayered  model  at  the  source  depth.  The 
dispersion  properties  that  control  the  far-field  waveform  are 
controlled  by  the  "global"  model  for  the  rest  of  the  path. 

The  points  made  in  the  previous  paragraphs  will  now  be 
clarified  by  considering  an  example.  As  an  average  path  model 
we  will  use  the  western  United  States  model  35-CM2  proposed  by 
Alexander  [1963],  which  is  tabulated  in  Table  1.  For  the 
source  region  we  use  the  same  model  modified  to  account  for 
details  of  the  local  geology.  Two  examples  will  be  considered. 
In  one  the  source  is  in  granite  like  PILEDRIVER  and  the  second 
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TABLE  1 

CRUSTAL  MODEL  35-CM2  FROM  ALEXANDER  [1963] 


Layer 


Depth 

(tan) 


Thickness 

(tan) 


a (km/sec)  B (tan/sec) 


1 

2.0 

2.0 

3.8 

1.75 

2.2 

2 

15.0 

13.0 

6.1 

3.6 

2.8 

3 

25.0 

10.0 

6.2 

3.7 

2.8 

4 

35.0 

10.0 

7.7 

4.1 

3.25 

5 

48.7 

13.7 

8.1 

4.65 

3.4 

6 

64.0 

15.3 

7.95 

4.55 

3.45 

7 

83.0 

19.0 

7.85 

4.5 

3.45 

8 

103.5 

20.5 

7.8 

4.4 

3.45 

CFUSTAL  MODEL  35-042  WITH  A PILEDRIVER  GRANITE  CAP 

i 1 

0.06 

0.06 

1.44 

0.83 

1.5 

( 2 

0.2 

0.14 

4.69 

2.71 

2.45 

3 

2.0 

1.8 

5.33 

2.78 

2.67 

4-9 

Layers 

2-8  of 

35-CM2 

CRUSTAL  MODEL  35-042  WITH  YUCCA  FLAT  TOFT  CAP 

1 

0.15 

0.15 

1.2 

0.8 

1.7 

2 

0.20 

0.05 

1.8 

1.0 

1.78 

3 

2.0 

1.8 

2.35 

1.3 

1.86 

4-9 

Layers 

2-8  of 

35-CM2 

source  area  is  Yucca  Flat.  The  two  source  structures  are  also 
tabulated  in  Table  1.  The  only  differences  between  the  three 
structures  are  in  the  top  two  kilometers. 

The  quantities  characterizing  mixed  path  models  made 
up  of  the  two  source  region  models  and  the  average  path  model 
of  Table  1 are  plotted  in  Figures  5,  6 and  7.  The  phase  and 
group  velocities  for  the  three  structures  are  shown  in  Figure 
5.  The  amplification  factor  AR  is  plotted  in  Figure  6.  In 


Figure  7 is  plotted  the  depth  excitation  factor,  K , for  the 
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two  source  region  models  and  the  transmission  coefficient,  T, 
for  the  boundaries.  From  Figures  5 and  6 we  see  that  the 
differing  material  in  the  top  two  kilometers  causes  substan- 
tial differences  in  the  quantities  plotted  for  periods  below 
10  to  15  seconds.  The  same  is  true  of  the  depth  excitation 
factor  in  Figure  7.  However,  the  transmission  coefficient 
is  within  20  percent  of  unity  for  periods  down  to  eight 
seconds . 


The  important  result  illustrated  by  the  foregoing 


examples  is  that  the  long  periods  at  which  M is  measured 

s 


pass  essentially  unchanged  through  material  boundaries  that 
are  restricted  to  the  top  few  kilometers  of  the  models. 

This  provides  a convincing  justification  for  computing  sur- 
face waves  according  to  the  formula 


A 
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aTslhA 
e / 
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where  r is  the  range  and  the  subscripted  quantities  are  as  in 
(8).  This  procedure  allows  us  to  use  a single  average  path 
model  for  events  in  a particular  source  region  while  account- 
ing for  changes  in  the  very  local  material  properties  in  a 
consistent  way. 
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AMPLIFICATION  FACTOR,  A (x  10u  microns/dyne) 


10  100 


PERIOD  (seconds) 


Figure  6.  The  amplification  factor  is  plotted  versus 
period  for  the  structure  of  Table  1. 


Figure  7.  The  transmission  coefficient  and  depth  excitation  factor  are  plotted 
versus  period  for  two  mixed  path  models. 
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4.3  THEORETICAL  SURFACE  WAVE  AMPLITUDE  - YIELD  RELATIONSHIP 

Using  the  theoretical  results  of  the  previous  section, 
theoretical  amplitude-yield  curves  are  constructed  for  explo- 
sions in  different  source  materials.  The  burial  depth  was 
assumed  to  vary  according  to  107  W^3  (meters) . Two  sets  of 
theoretical  seismograms  were  computed,  one  for  explosions  in 
granite  and  one  for  explosions  in  a material  similar  to  the 
saturated  tuffs  at  Yucca  Flat.  The  travel  path  models  were 
those  discussed  in  Section  4.2  as  examples  of  the  mixed  path 
technique.  They  are  tabulated  in  Table  1.  The  attenuation 
model  is  that  of  Tryggvason  [1965]  and  the  instrument  response 
is  for  the  standard  LRSM  long  period  system. 

The  theoretical  seismograms  are  shown  in  Figures  8 and 
9.  On  each  record  is  shown  the  depth  of  burial  and  yield  of 
the  source.  These  seismograms  are  for  a range  of  3000  km. 

For  the  periods  appearing  on  the  records  of  Figures  8 and  9 
the  source  function  ($ [w] ) is  essentially  constant  for  mate- 
rials of  interest.  Therefore,  the  only  important  characteristic 
of  the  source  potential  is  its  static  level  or  and  the 
seismograms  scale  with  this  quantity.  For  convenient  compari- 
son, all  theoretical  seismograms  of  this  section  were  computed 
with  the  same  source  function  which  is  characterized  by 
= 9.1  m3  at  a yield  of  0.02  kt. 

From  the  seismograms  of  Figures  8 and  9 we  see  that 
variations  in  yield  and  burial  depth  have  no  effect  on  the 
teleseismic  wave  form.  Comparing  the  two  sets,  it  is  clear 
that  changing  the  source  material  has  only  a very  minor 
effect.  Therefore,  it  is  easy  to  make  consistent  measure- 
ments for  computing  M_.  The  cycle  at  which  amplitude  mea- 

D 

surement  is  made  is  indicated  on  each  record  along  with  the 
Mb  which  is  computed  from 

I **“ 

Mg  — log  + log  A + 1.12  , (11) 


28 


MZonnHi  «zo»r>M3  nzonnwa  «2ojjom3  *>2oj iohj  moariMj  moafv-ij 


Yield 


Figure 


(km)  s 


. 2ft *00 


0 


-. 253*00 
0 

.745*00 


0 


-.437*00 

0 

.1*1*01 


0 


-.154*01 

0 

.304*01 


0 


-.25f*0i 

0 

.710*01 


0 


-.445*01 

0 

. 1f5*02 


0 


-.144*02 

0 

.334*02 


0 


.214*02 

0.  SO.  100.  ISO.  200.  250  . 300  . 350.  400.  4S0.  500.  550. 


10  0.23 

3.56 

A 

'wvvvVvuvl j 

lr  ■ 

50.  100.  ISO.  200.  250.  300.  350.  400.  450.  500.  550. 


. 25  0.31 

a A A A A A A AH 

, 3.95 

L 

50.  100.  150.  200.  250.  300.  350.  400.  450.  500.  550. 


100  0.50 

. u ^ ,/s  a A A A A A A A H 

A 4.56 

]\a 

^ vvvvvW\/\n 

Jlr  

50.  100.  150  . 200  . 250  . 300  . 350  . 400  . 450  . 500  . 550. 


250  0.67 

— r.  aaaAaAA 

AH 

llfl 4'.9.7  . 

■ ^ vVvVVVuv 

V | 

lr  — 

50.  100.  ISO.  200.  250.  300.  350.  400.  450.  500.  550 


1000  1.07 

^ ^ /\  A A A Ail  /I 

All 

ft,  5'S9 

- ^vvvVvvv 

tm 

lr 

. 

. Theoretical  seismograms  are  shown  for  a series  of 
explosions  in  NTS  granite.  The  yield,  burial 
depth  and  M_  are  shown  with  each  record  and  the 
cycle  at  which  Ms  is  measured  is  indicated  with  a 
bar.  The  period  of  this  cycle  is  19.7  seconds. 
Positive  vertical  is  down. 


A 


where  A is  peak-to-peak  amplitude  in  nanometers  and  A is  range 
in  degrees.  This  is  the  formula  used  by  Eisenhauer  [1976]. 

The  Mg  values  from  Figures  8 and  9 are  plotted  versus 
yield  in  Figure  10.  A linear  least  square  fit  was  made  to 
each  data  set  and  the  results  are  indicated  in  the  figure. 

The  first  thing  to  be  noticed  is  that  the  slope  is  slightly 
greater  than  unity.  Referring  to  (10)  , the  only  quantities 
that  change  from  one  yield  to  another  are  the  source  function 
(4')  , which  is  assumed  to  cube-root  scale  with  yield,  and  the 
depth  function  Kg^.  The  values  of  Kg^  at  20  seconds  are 
plotted  in  Figure  11  for  the  two  cases.  We  see  that  for  the 
granite  source  region  the  depth  function  is  nearly  constant 
over  this  range.  The  slope  of  1.02  is  then  due  to  the  source 
excitation  at  20  seconds  increasing  slightly  over  the  10  to 
1000  kt  yield  range.  The  greater  slope  of  1.06  for  the  tuff 
source  region  is  due  to  the  fact  that  the  Kg^  increases  some- 
what with  depth.  For  example,  with  all  other  factors  being 
equal,  an  explosion  in  tuff  at  a depth  of  1.07  km  will  have 
an  Mg  that  is  0.08  units  higher  than  one  at  a depth  of  0.23  km. 

The  second  interesting  feature  of  the  seismograms  of 
Figures  8 and  9 is  the  variation  of  amplitude  with  the  source 
material  while  the  source  function  (¥)  is  held  fixed.  From 
(10)  we  see  that  the  spectral  amplitude  ratio  is 


to 


(G) 
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where  the  superscripts  indicate  whether  the  source  is  in 
granite  (G)  or  tuff(T).  These  quantities  are  plotted  in 
Figures  5,  6,  and  7 or  tabulated  in  Table  1.  We  find  that 


Li. 


> 
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Then 


where  we  have  introduced  the  source  functions  which,  in  general 
are  important  for  the  scaling.  Referring  to  the  definition  of 
K (6) , we  see  that  y plays  an  important  role  in  this  param- 
eter.  In  view  of  this  fact,  we  find  it  useful  to  fit  the  ampli 
tude  ratio  by  the  relationship 


For  the  example  treated  here  the  n varies  from  0.83  at  10  kt 
to  0.71  at  1000  kt.  At  150  kt  we  have  n = 0.76.  A rough 
approximation  to  the  dependence  of  M on  the  source  and  source 

5 

material  properties  is  then 


where  ¥ is  the  (nearly  constant)  low  frequency  value  of  $ 
This  equation  should  be  approximately  be  valid  for  all  NTS 


source  regions  since  the  change  from  tuff  to  granite  in  the 
top  two  kilometers  of  the  path  model  (Table  1)  for  this 
example  is  about  the  most  severe  change  we  expect  to  encounter 


4.4  THE  DEPENDENCE  OF  M„  ON  EPICENTRAL  DISTANCE 

s 

The  Mg  values  presented  in  the  previous  section  were 
for  a particular  epicentral  distance  (3000  km) . In  that  sec- 
tion we  were  interested  in  the  change  in  M with  yield,  burial 

s 

depth  and  source  material.  In  subsequent  sections  the  absolute 
value  of  Mg  will  be  of  interest  and  we  should  determine  how 
well  this  quantity  is  defined  by  our  theoretical  analysis. 

In  Figure  12  we  show  seismograms  at  several  ranges  for 
the  case  in  which  the  yield  is  60  kt  and  the  depth  is  0.42  km 
in  granite.  For  each  of  these  records  the  amplitude  of  the 
cycle  with  period  closest  to  20  seconds  was  measured  and  Mg 
was  computed  according  to  two  formulas.  One  is  that  used  by 
Eisenhauer  [1976]  and  given  in  (11) . The  other  is  the  formula 
suggested  by  Marshall  and  Basham  [1972]  which  is 


Ms  “ log  I + B'(A)  + P<T>  ' 


(17) 


where  B'(A)  and  P(T)  are  tabulated  in  the  Marshall  and  Basham 
paper  (the  P(T)  appropriate  for  continental  North  America  was 
selected) . 

In  Table  2 the  Mfl  data  for  the  seismograms  of  Figure  12 
are  summarized.  The  Eisenhauer  formula  gives  Mg  values  that 
are  0.16  to  0.21  units  larger  than  those  from  the  other  formula. 
The  scatter  is  slightly  greater  for  the  values  from  the  Marshall 
and  Basham  formula.  The  3000  km  theoretical  record  gives  a 
theoretical  M that  is  close  to  the  "average"  value,  though  it 
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TABLE  2 


M VALUES  FROM  THE  SEISMOGRAMS  OF  FIGURE  11 
s 


Range 

(km) 

Amplitude 

(km) 

Period 

M 

s 

(Eisenhauer)  * 

P(T) 

M 

s 

(Marshall  & Basham) 

730 

4780 

18.5 

4.32 

-0.12 

4.14 

900 

4419 

18.9 

4.37 

-0.09 

4.21 

2000 

2123 

19.4 

4.40 

-0.05 

4.19 

3000 

1219 

19.7 

4.34 

-0.03 

4.13 

4000 

1085 

19.9 

4.41 

-0.01 

4.24 

Average  Mg 

4.37 

4.18 

We  have  corrected  for  the  instrument  response  at  the  actual 
period  rather  than  at  20  seconds,  which  may  or  may  not  be 
consistent  with  common  practice.  The  difference  this  may 
cause  is  small,  a few  hundredths  of  a unit  at  most. 
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is  perhaps  a bit  low.  It  seems  reasonable  to  use  the  M from 

s 

records  at  3000  km  as  our  theoretical  Mg  estimate. 

4.5  THE  DEPENDENCE  OF  ON  PROPERTIES  OF  THE  TRAVEL  PATH 

s 

For  computing  theoretical  surface  wave  seismograms  we 
are  using  two  crustal  models,  one  for  the  crust  in  the  vicinity 
of  the  source  and  one  for  the  "average"  crust  between  the 
source  and  receiver.  Even  if  the  source  were  azimuthally 
symmetric,  the  radiated  surface  waves  may  be  asymmetric  due 
to  azimuthal  differences  in  the  travel  path.  To  get  some  idea 
of  how  strong  this  effect  might  be,  we  computed  a series  of 
theoretical  seismograms  with  the  source  region  crust  being 
fixed  while  a number  of  different  continental  crustal  models 
were  used  for  the  rest  of  the  path. 

The  velocity-depth  profiles  for  the  various  crustal 
models  are  shown  in  Figure  13.  The  source  region  model  is 
one  proposed  by  McEvilly  [1964]  with  the  low  velocity  sur- 
face sedimentary  layer  replaced  by  PILEDRIVER  granite.  This 
model  is  much  like  the  Alexander  [1963]  model  35-CM2  with  a 
granite  cap  discussed  in  Section  3.2.  The  other  models  are: 


B - A North  American  crustal  model  from 
McEvilly  [1964] 

C - Model  35-CM2  from  Alexander  [1963] 

D - CIT  109  due  to  Archambeau,  et  al . , [1969] 
and  based  primarily  on  body  waves  from 
NTS  explosions 

E - A model  with  an  NTS  granite  crust  over  a 
mantle  derived  by  Helmberger  (private 
communication  from  Harkrider,  1977) 

F - A Gutenburg  continental  model  taken  from 
a paper  by  Harkrider,  at  al. , [1963] 

G - An  average  earth  model  constructed  by 

D.  Anderson  that  is  essentially  a monopole 
fit  to  the  free  oscillation  data. 
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The  phase  and  group  velocity  curves  for  these  structures  are 
shown  in  Figure  14 . 

As  is  apparent  from  the  formula  (10),  the  other  quantity 
of  importance  for  computing  surface  waves  is  the  transmission 
coefficient,  T(u).  Since  all  these  structures  are  pretty  much 
the  same,  this  coefficient  is  close  to  unity.  In  fact,  for 
models  B,  D,  E,  F and  G the  value  at  20  seconds  ranges  from  0.97 
to  1.03.  For  model  C,  the  average  earth  model,  the  structural 
contrast  is  greater  than  for  the  continental  models  and  T(w) 
at  20  seconds  is  0.87. 

Aside  from  the  T(w),  the  only  influence  of  the  average 
crust  model  (when  the  source  region  model  is  fixed)  is  in  the 
propagation  term,  the  Hankel  function  in  (10) . Then  when 
T(u)  is  nearly  unity  as  it  is  for  the  cases  discussed  here, 
the  amplitude  spectra  are  nearly  the  same  near  20  seconds  for 
all  the  path  models.  Differences  in  the  20  second  component 
in  the  time  domain  are  due  to  interference  effects. 

Synthetic  seismograms  at  3000  km  are  shown  in  Figure 
15  for  the  seven  path  models.  The  Mg  data  are  summarized  in 
Table  3 in  the  same  format  as  Table  2.  The  values  scatter 
by  a little  less  than  0.2  M.  units.  This  scatter  is  an  indica- 
tion  of  how  difficult  it  is  to  get  a good,  consistent  estimate 
of  the  amplitude  of  the  20  second  component  of  the  dispersed 
wave. 


4.6  THE  EFFECT  OF  SUPPRESSION  OF  THE  UPWARD  TRAVELING 
WAVES  GENERATED  BY  THE  SOURCE 

In  Section  3.2  we  outlined  several  potential  problems 
that  arise  when  the  explosion  is  represented  by  a spherically 
symmetric  source  of  elastic  waves  embedded  in  one  layer  of  a 
layered  earth  model.  We  noted  the  fact  that  for  most  events 
the  nonlinear  region  intersects  the  free  surface.  Also,  we 
know  the  source  region  is  inhomogeneous;  the  coupling 
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Figure  14.  The  Rayleigh  wave  phase  (left)  and  group  (right)  velocities  are  plotted 
for  the  structure  of  Figure  13. 


<*zo»nH3  Hioanx]  Mio^nni  «zo»nM3  MzoanMj  nzcam  >3  »zoxdh3 


STRUCTURE 


T 


.357*01 


-.337*01 

.327*01 


*.357*01 

.244*01 


-.240*01 
• 2?1*0t 


-.343*01 

.20**01 


-.m»oi 

.217*01 


-.214*01 

.1>1*01 


18.7 


19.9 


4.67 


4.55 


E 

--  a A A f' 

A A- 

u 

0.9  4.65 

733 

^ — 

VIM 

IT 

i/vWf' 

jyi/wvw1*  ■ , 

. . i 

4.70 


-.1***01 


Figure  15. 


Synthetic  seismograms  are  shown  for  a 100  kt  explo- 
sion at  a depth  of  0.7  km,  using  a mixed  path  model. 
The  source  region  model  is  A and  seven  average  path 
models  are  used.  The  time  corresponding  to  zero  on 
the  scale  is  indicated  on  each  record  as  well  as  the 
M?  (using  (11))  and  Ms  period.  A small  arrow  on  the 
time  axis  marks  900  seconds  after  detonation  time. 
The  cycle  at  which  the  Mr  measurement  is  made  is 
marked  with  a bar.  At  the  left  is  amplitude  in 
microns  at  25  seconds. 
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properties  of  materials  above  and  below  the  source  are  not  the 
same.  The  most  obvious  difference  is  the  fact  that  the 
"average"  overburden  pressure  in  the  inelastic  region  above 
the  source  is  lower  than  in  the  analogous  region  below  the 
source.  Lower  overburden  pressures  generally  lead  to  larger 
source  functions.  Acting  in  the  opposite  direction  is  the 
effect  of  the  generally  higher  air-filled  porosity  above  the 
source.  Perhaps  more  important  is  the  fact  that  at  least 
some  energy  is  lost  from  the  upgoing  wave  due  to  spallation. 
The  upgoing  wave  would  also  seem  to  be  more  subject  to  losses 
due  to  scattering  and  inelastic  attenuation  in  the  immediate 
vicinity  of  the  source  than  the  downgoing  wave.  At  present 
we  do  not  really  know  how  important  these  various  effects  are. 
We  need  to  carry  out  some  theoretical  studies  involving  two- 
dimensional  finite  difference  calculations  to  develop  a better 
understanding.  Even  so,  with  our  relatively  simple  models  we 
can  do  some  things  to  better  delineate  the  effects  these 
source  asymmetries  might  have. 

Our  theoretical  formulation  for  both  body  and  surface 
waves  has  been  modified  to  allow  total  or  partial  suppression 
of  the  source  radiation  emitted  upward  or  downward  from  the 
source  (Harkrider  and  Bache  [1977]).  This  allows  us  to  do 
several  things  to  better  understand  the  effects  being  dis- 
cussed. First,  we  can  hypothesize  that  the  waves  traveling 
upward  from  the  source  (takeoff  angles  above  the  horizontal) 
are  trapped  or  scattered  or  in  some  way  prevented  from 
coupling  into  the  far-field  Rayleigh  wave.  We  can  then  com- 
pare the  Rayleigh  wave  amplitudes  or  M from  the  total  source 
to  those  from  the  same  source  with  the  upward  traveling  waves 
suppressed.  A second  useful  thing  to  do  is  to  suppose  that 
the  equivalent  elastic  source  (¥  (w) ) is  different  for  the  up- 
going  and  downgoing  waves.  In  fact,  since  the  geologic 
environment  is  generally  different  above  and  below  the  source. 
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we  can  obtain  first  order  estimates  of  how  different  these 
sources  might  be  with  relatively  inexpensive  one-dimensional 
source  calculations.  Even  without  such  source  calculations, 
we  can  delineate  the  range  of  possible  effects  since  we 
believe  the  source  function  is  nearly  constant  over  the 
frequency  band  of  interest  for  surface  waves.  Thus,  one- 
dimensional sources  can  be  represented  by  a single  number, 

Y*.  A simple  way  to  investigate  the  question  is  then  to 
suppose  that  the  size  of  the  appropriate  for  the  upgoing 
waves  is  different  than  that  for  the  downgoing  waves. 

Let  us  return  to  the  calculations  of  Section  3.3  for 
which  theoretical  seismograms  are  shown  in  Figures  8 and  9. 

The  question  is,  what  is  the  effect  of  totally  or  partially 
suppressing  the  upgoing  waves  excited  by  the  source?  This 
effect  is  delineated  for  one  case  by  the  seismograms  shown 
in  Figure  16.  The  example  studied  is  the  granite  source, 

60  kt  yield,  0.42  km  depth  case  from  Figure  8.  We  show  the 
seismograms  for  the  upgoing  and  downgoing  waves  separately 
and  for  their  total,  which  is  the  same  as  the  corresponding 
record  in  Figure  8.  Seismograms  computed  with  the  source 
{V'n)  for  the  upgoing  waves  being  0.25,  0.50  and  0.75  times 
the  size  of  the  source  for  the  downgoing  waves  are  also 
shown.  For  this  case  the  surface  waves  for  the  up  and  down- 
going waves  are  almost  exactly  in  phase.  The  addition  of 
various  fractions  of  the  upgoing  wave  portion  of  the  source 
then  has  rather  simple  and  predictable  effects. 

In  Figure  17  we  study  the  contribution  of  the  up  and 
downgoing  portions  of  the  source  for  an  event  just  like  that 
in  Figure  16  except  that  the  source  is  in  tuff.  For  this  case 
the  seismogram  for  the  upgoing  portion  of  the  source  is 
nearly  the  negative  of  that  for  the  downgoing  portion  with 
a slight  shift  to  later  time.  Then  in  contrast  to  the  resuics 
for  the  source  in  granite,  addition  of  the  two  portions  of  the 

source  leads  to  destructive  interference  and  a reduction  in  M . 
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Figure  16. 


Seismograms  are  shown  for  the  60  kt  granite  source 
at  0.42  km  depth  from  Figure  8.  The  upgoing  wave 
and  downgoing  wave  contributions  to  the  source  are 
shown  separately.  The  seismograms  for  the  down- 
going wave  portion  added  to  the  indicated  fraction 
of  the  upgoing  portion  are  also  shown.  Note  that 
the  polarity  is  reversed  from  earlier  seismograms 
(positive  vertical  is  up) . 


Figure  17.  Seismograms  are  shown  that  are  identical  to  those 
in  Figure  16  except  that  the  source  is  in  tuff. 
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While  some  part  of  the  Mg  variation  in  Figure  17  is  due  to 
shifts  in  the  phase  at  which  the  amplitude  is  measured,  com- 
parison of  the  20  second  spectral  amplitudes  for  these  seismo- 
grams gives  much  the  same  result. 

What  is  the  effect  of  suppressing  the  upward  traveling 
waves  on  the  Mg-yield  curves?  In  Figures  18  and  19  we  show 
seismograms  for  the  same  events  as  in  Figures  8 and  9 except 
that  the  upward  traveling  waves  are  entirely  suppressed.  As 
was  the  case  in  Figures  8 and  9,  we  note  that  changes  in  burial 
depth  over  this  range  have  essentially  no  effect  on  the  wave- 
form. The  M values  are  plotted  versus  yield  in  Figure  20 
together  with  least  square  fits  to  these  data  and  the  analogous 
data  from  Figure  10.  The  only  really  significant  effect  of 
suppressing  the  upgoing  waves  is  that  noted  in  connection  with 
Figures  16  and  17;  the  M is  reduced  for  the  granite  source  and 
is  enhanced  for  the  tuff  source. 

Let  us  now  summarize  what  we  have  learned  from  the 
examples  of  this  section.  It  seems  quite  plausible  to  suppose 
that  the  upgoing  wave  is  suppressed  to  some  degree  compared 
to  the  downgoing  wave  for  many  explosions.  For  a source  in 
a low  strength  material  overlying  a hard  basement,  we  find 
that  the  surface  wave  excited  by  the  upgoing  waves  is  about 
225°  out  of  phase  with  that  associated  with  the  downgoing 
wave.  Suppression  of  the  upgoing  wave  then  acts  to  increase 
M compared  to  that  for  the  total  source.  On  the  other  hand, 

S 

when  the  source  is  in  hard  rock,  the  two  portions  of  the  sur- 
face wave  are  in  phase  and  suppression  of  the  upgoing  wave 
reduces  the  Mg.  Recall  that  we  found  (Section  4.3,  Equation 

(16))  for  the  total  source  that  M a log  1 . When  the 

s s 

upgoing  wave  is  totally  or  partially  suppressed,  the  depen- 
ence  on  yg  is  much  weaker.  In  fact,  for  total  suppression  of 
the  upgoing  wave  we  find  that 

Mg  a log  [pg/8  . (18) 
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Figure  18.  Theoretical  seismograms  are  shown  that  are 
identical  to  those  of  Figure  8 (the  source 
is  in  granite)  except  that  the  upgoing 
waves  are  suppressed.  Positive  vertical 
is  down. 


.*<2-01 


Yield 

10 


Dep 
(km) 

0.23 


s 

3.07 


.755  01  J — 1 1 

0.  50.  100.  ISO.  200.  2S0.  000.  350.  <00.  <50.  500.  550. 

213,00  r j 

25  0.31  ft  3.47 


. Uf.00  * * — 1 — 1 

0.  50.  100.  ISO.  200.  250.  300.  350.  <00.  <50.  500.  550. 
.523*00  1 

60  0.42  ft  3.85 


. <57*00  * * * — 1 * * * 

0.  50.  100.  150.  200.  250.  300.  350.  <00.  <50.  500.  550. 

,**s*oo  i 1 1 


100  0.50 

^ /-\  /\  A A A A 

/L 

J 

4.07 

^ V V V V \j 

TT 

\/V 

V 

.745*00  * * * * — — * * * 

0.  SO.  100.  150.  200.  250.  300.  350.  <00.  <50.  500.  550. 
.230*01  j 

250  0.67  ft  4.48 


0.  50.  100.  150.  200.  250.  300.  350.  <00.  <50.  500.  550. 


0.  50.  100.  ISO.  200.  250  . 300.  350.  <00.  <50.  500.  550. 

, 1 02*02  -j 

1000  1.07  ft  5.10 


0 . 50.  100.  150.  200.  250  . 300  . 350.  <00.  <50.  500.  550. 


Figure  19. 


Theoretical  seismograms  are  shown  that  are 
identical  to  those  of  Figure  9 (the  source 
is  in  tuff)  except  that  the  upgoing  waves 
are  suppressed.  Positive  vertical  is  down. 
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That  is,  in  this  case  the  Mg  is  nearly  independent  of  the 
properties  of  the  source  medium  other  than  indirectly  via 

nw) . 


V.  SIGNAL  DECONVOLUTION 


5 . 1 INTRODUCTION 

Teleseismic  signals  are  controlled  by  characteristics 
of  the  source  and  the  travel  path.  Since  most  underground 
nuclear  tests  are  conducted  in  a small  number  of  known  test 
sites,  there  are  fairly  large  populations  of  events  for  which 
the  path  effects  should  be  about  the  same.  Even  comparing 
events  at  different  test  sites,  local  station  effects  might 
be  expected  to  be  stationary.  In  short,  it  should  be  possible 
to  isolate  some  components  of  the  transfer  function  relating 
the  explosion  source  to  the  teleseismic  waveform  and,  thereby, 
improve  our  ability  to  isolate  the  source  effects. 

To  be  more  explicit,  let  us  express  the  teleseismic 
short  period  signal  in  the  following  form: 

u - TI  • trc  • tm  • Tsc  • tes  • M ' (19) 

where  W is  the  yield  of  the  explosion  and  the  T represent  the 
transfer  functions  that  relate  the  explosion  to  the  recorded 
waveform.  The  TEg  represents  the  conversion  of  the  energy  of 
the  explosion  into  elastic  waves  in  the  earth  and  depends  on 
the  local  materal  properties.  The  most  one  could  hope  for 
with  seismic  methods  is  to  estimate  the  equivalent  elastic 
source,  or  TEg  • W.  The  Tgc,  TM  and  TRC  break  the  travel 
path  into  three  segments,  the  crust  at  the  source,  the  upper 
mantle  and  the  crust  at  the  receiver.  The  Tj  is  the  transfer 
function  of  the  seismometer  and  is  easy  to  remove. 

How  can  we  "deconvolve"  recorded  seismograms  to  isolate 
and  quantify  these  transfer  functions?  There  are  several 
schemes  that  come  to  mind: 

1.  Consider  a series  of  teleseismic  recordings 
of  the  same  event.  The  common  part  should 
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be  Tj  • Tgc  • Tes  • W.  If  this  can  be 


extracted  and  deconvolved  with  the  seis- 
mograms, the  residual  should  be  TM  • TRC. 

We  will  know  how  well  this  works  by  how 
constant  TM  • TRC  is  from  event-to-event . 

2.  Consider  a series  of  similar  events  recorded 
at  a particular  station.  The  common  part 
should  be  • TRC  • TM  and  the  residual  is 
then  Tgc  • TRS  • W.  The  Tj  • TRC  • T^  should 
be  nearly  constant  from  event-to-event. 

3.  Having  done  previous  study,  we  can  further 
analyze  the  T^.  • TRC  • T^  to  determine  TRC< 

The  Tj  is  known  and  the  TM  should  be  nearly 
independent  of  frequency  for  most  stations 
beyond  35°  or  40°.  We  then  can  estimate  TRC 
for  the  particular  station  under  study.  We 
hope  this  TRC  would  be  fairly  constant  from 
event-to-event,  though  it  may  exhibit  some 
range  (incidence  angle)  and  azimuthal 
dependence. 

We  have  outlined  what  we  would  like  to  be  able  to 
achieve  by  deconvolving  teleseismic  short  period  recordings 
of  underground  explosions.  In  the  remainder  of  this  section 
we  outline  the  method  we  plan  to  use  which  is  closely  related 
to  the  homomorphic  deconvolution  method  described  by  Clayton 
and  Wiggins  [1976]. 

Many  schemes  have  been  suggested  for  inverting 
seismic  data  to  isolate  the  signal  that  is  common  to  all  the 
observations.  Problems  of  this  nature  belong  to  the  type 
designated  through  various  names:  deconvolution,  inverse 
filtering,  decomposition,  etc.  All  these  operations,  being 
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of  the  inverse  type,  are  characterized  by  their  potential 
instabilities.  In  this  respect,  they  are  all  similar  to 
the  simple  inverse  arithmetical  operation,  i.e.,  the  divi- 
sion, the  inverse  of  the  multiplication  operation.  If  the 
divisor  gets  to  be  very  small,  the  inverse  grows  to  very 
large  size  and  greatly  magnifies  the  quantity  on  which  it 
operates.  If  this  quantity  is  not  known  with  sufficient 
accuracy  or  is  contaminated  with  noise,  the  result  of  the 
operation  becomes  meaningless. 

Deconvolution  or  decomposition  of  a given  signal  is 
entirely  controlled  by  the  assumptions  that  the  user  of  the 
technique  is  willing  to  make.  If  these  assumptions  are  valid 
and  the  data  fits  the  model  they  imply,  deconvolution  can  be 
a powerful  tool  indeed;  if,  on  the  contrary,  they  are  irrele- 
vant, the  results  may  be  very  discouraging  or  even  disastrous. 

Estimating  a source  signature  is  equivalent  to  per- 
forming the  decomposition  of  an  observation,  sometimes  called 
an  output,  into  its  two  basic  components:  a waveshape,  also 
called  the  input,  and  a response  or  transfer  function. 

Response  or  transfer  function  analysis  originated 
primarily  in  the  electrical  engineering  field  and,  in  that 
field,  it  has  branched  out  in  several  separate  expertises 
such  as  network  analysis,  filter  design,  time  series  analysis, 
information  and  communication  theory,  control  theory,  etc. 

More  recently,  it  has  been  extensively  applied  and  adapted  to 
the  seismic  exploration  area,  particularly  with  the  intro- 
duction of  optimal  filtering  and  deconvolution  techniques 
in  the  early  fifties. 

In  exploration  reflection  seismology,  most  deconvolu- 
tion techniques  rely  on  the  assumption  that  there  are  enough 
events  of  constant  shape  in  a single  trace  over  the  time 
interval  of  interest  that  the  theory  of  stationary  time  series 
applies.  This  is  equivalent  to  saying  that  the  transfer 
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function  or  response  is  white  while  the  wave  shape  has 

: 

minimum  phase  characteristics.  With  these  two  assumptions, 
the  decomposition  is  unique  and  the  answer  can  be  obtained 
from  a single  autocorrelation  analysis.  The  autocorrelation 
of  the  signal  determines  the  power  spectrum  from  which  the 
minimum  phase  waveshape  can  be  derived  and  its  inverse  used 
to  perform  the  deconvolution  of  the  trace.  Unfortunately, 
these  assumptions  appear  unrealistic  in  the  case  of  tele- 
seismic  events.  We  must  then  resort  to  a different  approach. 

Another  approach  to  deconvolution  was  developed  in 
the  late  sixties  by  generalizing  the  notion  of  superposition 
through  the  introduction  of  a nonlinear  homomorphic  trans- 
formation, the  logarithmic  operation.  Just  remember  how, 
through  the  logarithmic  operation,  the  arithmetical  operation 
of  multiplication  is  reduced  to  an  addition,  a very  substan- 
tial simplification.  Similarly,  when  considering  the  con- 
volution operation,  note  that  in  the  frequency  domain  it 
reduces  to  a complex  multiplication.  Therefore,  if  one 
makes  a logarithmic  transformation  in  the  spectral  domain, 
the  multiplication  reduces  to  a mere  addition.  This  sounds 
too  good  to  be  true  and,  in  fact,  we  are  going  to  show  that 
it  is  not  quite  that  simple,  but  it  is  worth  pursuing. 

The  spectrum  of  a signal  is  generally  a complex 
quantity  and  the  logarithm  of  a complex  quantity  is  not 
single  valued,  since  if  z = pe^6  = pe^0+^KlT^  An  z = in  p + 
i(0+2Kir)  where  K is  any  positive  or  negative  integer.  Thus, 
we  have  an  inherent  ambiguity  in  using  the  spectral  route. 

However,  in  a way  similar  to  the  minimum  phase  assumption  in 
the  ordinary  deconvolution,  we  may  be  able  to  remove  the  2n 
ambiguity  by  making  again  further  assumptions. 

This  deconvolution  technique,  called  homomorphic 
deconvolution,  has  been  used  by  several  teleseismic  analysts 
who  think  that  it  is  better  adapted  to  the  processing  of 
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teleseismic  recordings  than  the  deconvolution  techniques 
employed  in  seismic  exploration.  They  claim  it  is  partic- 
ularly well  adapted  to  echo  removal  or  dereverberation. 

In  order  to  see  why,  we  must  first  describe  more  precisely 
what  homomorphic  deconvolution  is. 

Assume  a given  signal,  s (t) , it  is  known  to  be  the 
result  of  the  convolution  of  two  other  signals,  a waveshape 
w(t)  and  a response  r(t).  Then 

t t 

s(t)  = w(t)*r(t)  = f w(t)r(T-t)  dx  = [ r(t)w(x-t)  dx  . (20) 


In  the  frequency  domain  this  translates  into 


S ( oo ) = W(oj)  • R(u)  = | S (to)  | el6(u)  = |w(u)  | • |R(u)|  el(4>(a,)  + (6)) 


(21) 


where  capital  letters  are  used  to  represent  Fourier  transforms 
of  the  corresponding  signal  which  is  denoted  by  a small  letter. 
Taking  the  logarithm,  we  have 


£n  S (u)  = £n  |S(u>)|  + 16  (w)  = £n  |W(u)  | + £n|R(w) 


+ i [ 4>  (u)>  + Y(u)  + 2Ktt]  . 


(22) 


A new  quantity  is  then  introduced  by  taking  the  inverse 
Fourier  transform  of  the  £n  S(u).  This  produces  a signal  which, 
while  it  is  real,  has  been  named  the  complex  cepstrum  to  differ- 
entiate it  from  the  real  cepstrum  (e.g.,  reference)  which  uses 
the  same  transformation  on  the  power  spectrum,  a real  quantity. 
The  real  cepstrum  is  in  the  time  domain,  but  to  retain  the  dis- 
tinction between  it  and  the  original  input  signal,  this  new 
domain  is  called  the  quefrency  domain. 
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Our  technique  for  a signal  deconvolution  takes  the 
following  form.  Given  a number  of  signals  having  a good 
degree  of  similarity,  we  first  take  their  Fourier  transform 
to  derive  for  each  signal  its  amplitude  spectrum  and  its 
phase  in  the  usual  (-it, it)  interval.  The  phase  is  then  un- 
wrapped in  a very  simple  way.  If  the  phase  moves  from  the 
second  to  the  third  quadrant,  2tt  is  subtracted.  If  it  goes 
from  the  third  quadrant  to  the  second,  2tt  is  added.  In 
this  way,  we  obtain  a reasonably  smooth  phase  if  the  sampling 
rate  produces  a high  enough  Nyquist  frequency.  Actually,  the 
criterion  used  is  whether  * $(i+l)  - $(i)  is  greater  than 
3ir/2  or  smaller  than  -3ir/2.  When  the  phase  moves  from  the 
second  quadrant  to  the  third,  -tt  < $ (I)  < - ir/2  and 
ir/2  <_  $(I+1)  <,  so  that  tt  <_  <_  2ir.  Similarly,  when  the 

phase  moves  from  the  third  quadrant  to  the  second,  ir/2  <_  ^ (I)  £ tt 
and  -it  <_  $(I+1)  <_  -ir/2,  so  that  -2tt  <_  £ -it.  The  selection 

of  3ir/2  equally  distributes  the  probabilities.  Another 
approach  would  be  to  add  and  subtract  2tt  to  the  value  and 
select  the  one  of  the  three  choices  which  has  the  smallest 
absolute  amplitude.  After  the  phase  has  been  unwrapped,  the 
logarithmic  transformation  is  performed  to  produce  a complex 
quantity  (in  S(u))  for  each  signal. 

The  next  step  is  to  arithmetically  average  the  amplitude 
and  phase  of  the  logarithm  of  the  transform  for  each  signal. 

That  is,  if  we  had  N signals  with  transforms 

. , , . i>|> . (u» 

S.(u>)  * |w(u)  | |R±  (co)  | e1*  w e 1 , i - 1,  N,  (23) 

then 

in  S^(u)  = in|w(w)  | + in|R^(w)|  + i [4>(w)  + ^(w) 

(24) 


+ 2Kir  + 2J.Tr)  , 


j 


where  K and  J\  are  integers.  The  averaging  is  then  done 
according  to 


A(u>) 


Y Un|W(«)|  + in |R. (u) | J 
5=1 


£ [<(>(10)  + ^i(u)  + 2Ktt  + 2Jiir] 


B(o>)  = — 


with  the  average  logarithmic  spectrum  being 


X (u)  = A (to)  + iB  (u)  . 


The  spectrum  of  the  average  source  shape  is  then 


Y (w)  - eA(“>  eiB<u) 


The  agreement  between  Y(w)  and  the  common  wave  shape,  W(w),  then 
depends  on  how  symmetrically  the  (w)  • W(u)  are  distributed 
about  the  W(w). 

For  a final  estimate  of  the  average  signal,  it  is  often 
advantageous  to  truncate  or  smooth  the  spectrum  Y(u).  We  are 
currently  applying  a cosine  taper  to  the  high  frequency  portion 
of  Y(w).  Our  results  indicate  that  most  of  the  important 
information  is  contained  in  the  low  frequency  portion  of  the 
spectrum  and  rather  severe  low  pass  filters  can  be  applied 
without  seriously  affecting  the  solution. 

Estimating  a signal  shape  is  tricky  business  requiring 
adequate  understanding  of  the  various  factors  controlling  the 
presence  of  the  signal  shape  in  the  individual  data  inputs. 
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A common  problem  is  that  of  time  relationship,  another  is 
that  of  noise  contamination  and  still  a third  one  is  com- 
plexity of  the  source.  An  adequate  procedure  must  make 
optimum  use  of  the  information  contributed  by  each  input 
signal.  The  scheme  we  have  just  described  appears  to  do 
that  adequately,  probably  because  the  spectral  averaging 
through  the  logarithmic  transformation  leads  to  a geometric 
average  instead  of  an  arithmetic  one,  thus  reducing  even 
more  the  variations  from  the  average. 

In  Section  5.1  we  set  forth  the  criterion  by  which 
the  adequacy  of  this  technique  is  to  be  judged;  that  is, 
by  the  degree  by  which  the  transfer  functions  remain  the 
same  from  one  set  of  data  to  another.  To  make  this  judg- 
ment, we  need  to  test  the  technique  on  a representative 
sampling  of  data  and  such  tests  are  now  underway.  As  an 
indication  of  the  results,  a few  simple  examples  are  shown 
in  the  following  section. 

5.2  PRELIMINARY  RESUL1C 

The  first  sampling  of  seismograms  analyzed  consists 
of  four  recordings  of  an  intermediate  yield  explosion  at 
Pahute  Mesa,  NTS . The  event  is  MAST  and  the  recordings  are 
shown  in  Figure  21.  At  the  bottom  of  the  figure  are  two 
examples  of  the  average  wavelet  present  in  these  data. 

First  shown  is  the  average  of  all  four  seismograms.  Then 
station  EBC  was  dropped  and  the  remaining  three  traces  were 
averaged.  Comparing  the  two,  we  get  some  indication  of  how 
much  the  average  time  signal  can  vary  by  increasing  or 
decreasing  the  number  of  traces  used  in  the  average.  Of 
course,  the  more  pertinent  question  is  how  sensitive  the 
transfer  function  at  individual  stations  is  to  this  average. 

A second  example  is  shown  in  Figure  22.  In  this  case 
we  have  seven  recordings  of  a presumed  Soviet  explosion.  The 
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JU01 


ESPZH 


EIM 


4 STATION 
AVERAGE 


3 STATION 
AVERAGE 


Figure  21.  Comparison  of  average  wavelets  determined  from 
MAST  seismograms  at  4 Stations  (EBC,  JU01, 
ESPZH  and  EIM)  and  3 Stations  (JU01,  ESPZH  and 
EIM)  . 
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interesting  thing  about  the  average  signal  is  that  it  is  a 
quite  discrete  pulse  having  about  the  shape  we  expect  for  a 
buried  explosion.  Once  again,  we  must  say  that  we  need  to 
compute  more  examples  and,  especially,  study  the  transfer 
functions,  to  quantify  how  good  this  average  pulse  is. 

In  summary,  we  have  developed  a signal  averaging  or 


decomposition  scheme  that  is  rather  simple  to  implement  and 
apply.  The  crucial  operation  in  applying  this  technique  is 
the  unwrapping  of  the  phase  of  the  spectrum  of  the  individual 
signals.  A few  very  preliminary  examples  show  that  the  tech- 
nique derives  an  average  signal  that  seems  to  be  reasonable 
for  the  data  being  processed.  The  next  step,  which  is  now 
underway,  is  to  compute  transfer  functions  from  which  the 
validity  and  utility  of  the  technique  can  be  judged  according 
to  the  criteria  set  forth  in  the  introduction  to  this  section 
(Section  5.1). 
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